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ABSTRACT 

We apply the Om diagnostic to models for dark energy based on scalar fields. In case of 
the power law potentials, we demonstrate the possibility of slowing down the expansion 
of the Universe around the present epoch for a specific range in the parameter space. 
For these models, we also examine the issues concerning the age of Universe. We 
use the Om diagnostic to distinguish the ACDM model from non minimally coupled 
scalar field, phantom field and generic quintessence models. Our study shows that the 
Om has zero, positive and negative curvatures for ACDM, phantom and quintessence 
models respectively. We use an integrated data base (SN+Hubble+BAO+CMB) for 
observational analysis and demonstrate that Om is a useful diagnostic to apply to 
observational data. 

Key words: cosmological parameters— cosmology: observations— cosmology: 
theory— dark energy. 


1 INTRODUCTION 

A number of observational investigations such as, Type la 
supernovae (Riess et al. 1998; Perlmutter et al. 1999), cos¬ 
mic microwave background radiation (Spergel et al. 2003; 
Komatsu et al. 2011), surveys of the large scale structure 
(Eisenstein et al. 2005) and PLANCK 2013 (Ade et al. 2013) 
indicate that our Universe is accelerating at present. In the 
standard framework based upon Einstein gravity, cosmic ac¬ 
celeration can be explained by an exotic fluid with large 
negative pressure filling the Universe, dubbed 1 dark energy’ 
(Riess et al. 1998; Perlmutter et al. 1999; Sahni & Starobin- 
sky 2000; Spergel et al. 2003; Copeland, Sami & Tsujikawa 
2006; Sami 2009; Komatsu et al. 2011; Ade et al. 2013; Sami 
& Myrzakulov 2013). The simplest candidate for dark en¬ 
ergy (DE)— the cosmological constant A is plagued with 
difficult theoretical issues related to fine tuning. A variety of 
other dark energy models have been studied in the literature, 
namely, quintessence (Ratra & Peebels 1988), ACDM (Wein¬ 
berg 1989) and the phantom field (Caldwell, Kamionkowski 
& Weinberg 2003; Setare 2007). 

Alternatively, cosmic acceleration can also be mimicked 
by the large scale modification of gravity (Boisseau et al. 
2000; Dvali, Gabadadze & Porrati 2000; Sahni & Shtanov 
2003; Sahni, Shtanov & Viznyuk 2005; Trodden 2007; Bre- 
vik 2008; Ali, Gannouji & Sami 2010; Gannouji & Sami 
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2010). Though, the late time cosmic acceleration is consid¬ 
ered to be an established phenomenon at present, its under¬ 
lying cause remains uncertain. A large number of models, 
within the framework of standard lore and modified theo¬ 
ries of gravity can explain the said phenomenon. It is there¬ 
fore important to find ways to discriminate between var¬ 
ious competing models. To this effect, important geomet¬ 
rical diagnostics have been recently suggested in the liter¬ 
ature such as, Statefinder (Alam et al. 2003; Sahni et al. 
2003), Statefinder hierarchy (Arabsalmani & Sahni 2011), 
Om (Sahni, Shafieloo & Starobinsky 2008, 2014) and Om3 
(Shafieloo, Sahni & Starobinsky 2012). Statefinders use the 
second, third and higher order derivatives of the scale factor 
with respect to cosmic time whereas Om relies on first order 
derivative alone. Consequently Om is a simpler diagnostic 
when applied to observations. 

The statefinder method has been extensively used in the 
literature to distinguish among various models of dark en¬ 
ergy and modified theories of gravity. For instance, non min¬ 
imally coupled scalar field, galileon field, Dvali, Gabadadze 
and Porrati (DGP) model, bimetric (bigravity) theory of 
massive gravity and others have been investigated (Sami 
et al. 2012; Myrzakulov & Shahalam 2013) using this di¬ 
agnostic. Observational constraints have been put on the 
statefinder pair {r, s} and deceleration parameter q using 
Union2.1 compilation data (Suzuki et al. 2012) and 28 points 
of Hubble data (Farooq & Ratra 2013), for the power law 
cosmological model (Rani et al. 2014). The statefinder anal¬ 
ysis has been applied to non minimally coupled galileons 
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and f(T) cosmology by several authors (Jamil et al. 2012; 
Jamil, Momeni & Myrzakulov 2013). Recently, statefinder 
hierarchy has been applied to models based upon Chaplygin 
gas, light mass galileons and holographic dark energy (Li, 
Yang & Chen 2014; Myrzakulov & Shahalam 2014; Zhang, 
Cui & Zhang 2014). 

Coming back to Om diagnostic, to be employed in this 
paper, we should note its excellent features. First, it can 
discriminate dynamical dark energy models from ACDM, 
in a robust way, even if the value of the matter density is 
not precisely known. Secondly, it can provide a null test of 
ACDM hypothesis, i.e., Om(z) - f2o m =0, if dark energy is a 
cosmological constant. Om has zero, negative and positive 
curvatures for ACDM, quintessence and phantom models 
respectively. 

Our paper is organised as follows: In section[2j we briefly 
revisit the Om diagnostic, to be used in the subsequent sec¬ 
tions. In section [21 we display evolution equations of phan¬ 
tom (non-phantom) scalar field in the autonomous form. In 
sub-section l3.ll we consider power law potentials and exam¬ 
ine the slowing down of cosmic acceleration for a possible 
range of the parametric space, whereas in sub-section l3.2l we 
consider age of the Universe in these models. In sub-section 
E3 we examine equation of state (EOS) and Om behaviour 
for the tracking potential V(tj>) = Vo [cosh (a(j>/M p ) — l] p . 
Sub-section and section [4] are devoted to application of 
Om to phantom field with linear potential and non mini¬ 
mally coupled scalar field respectively. Our results are sum¬ 
marized in the last section. In appendix A, we carry out joint 
data analysis and put observational constraints on the model 
parameters using Union2.1 compilation data (Suzuki et al. 
2012), 28 points of Hubble data (Farooq & Ratra 2013), 
BAO data (Eisenstein et al. 2005; Percival et al. 2010; Beut- 
ler et al. 2011; Blake et al. 2011; Jarosik et al. 2011; Giostri 
et al. 2012) and CMB data (Komatsu et al. 2011). 


2 Om DIAGNOSTIC 


Om is a geometrical diagnostic which combines Hubble pa¬ 
rameter and redshift. It can differentiate a dynamical dark 
energy model from ACDM, with and without reference to 
matter density. Constant behaviour of Om{z) with respect 
to 0 signifies that DE is a cosmological constant (ACDM). 
The positive slope of Om(z) implies that dark energy is 
phantom (w < — 1) whereas the negative slope means that 
DE behaves like quintessence (w > —1). Following Sahni et 
al. (2008); Zunckel & Clarkson (2008), Om(z) for spatially 
flat Universe is defined as 


Om(z) 


H 2 (z)/H 2 -1 
(l + 2) 3 -l 


( 1 ) 


where Ho is the present value of the Hubble parameter. 
Since, for constant equation of state parameter, 

H 2 {z) = jy 0 2 (fl 0 m(l + z ) 3 + (1 - n 0 m)(l + z) 3(1+W) ), (2) 


we have the following expression for Om(z), 

(l 4. _\ 3 (i+w) _ i 

Om(z) = n 0 m + (1 - n 0m ) y ; —— , (3) 

(1 + z)* - 1 
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Figure 1. This figure shows the evolution of Om(z) versus the 
redshift z for DE models with flom = 0.3 and w = -1.5, -1.1, 
-1, -0.7, -0.4, -0.1 (lower to top). The horizontal line represents 
ACDM with w = -1, and has zero curvature. The DE models 
with w > — 1 (quintessence) have negative curvature whereas DE 
models with w < — 1 (phantom) have positive curvature. 


for phantom (w < —1). This kind of behaviour can be clearly 
seen in figureQ] It is remarkable that Om provides a null test 
of the ACDM hypothesis. Another conclusion, that follows 
from this figure is that the growth of Om(z) at late time 
favours the decaying dark energy models (Shafieloo, Sahni 
& Starobinsky 2009). 


3 SCALAR FIELD DYNAMICS 

In this section, we will examine cosmological dynamics of 
phantom (non-phantom) field with generic potentials, using 
Om diagnostic. We will be interested in distinguishing fea¬ 
tures of the dynamics which allow us to differentiate these 
models from ACDM. Secondly, we will investigate the possi¬ 
bility of slowing down of cosmic acceleration in these mod¬ 
els. The models showing transient acceleration have been 
studied previously (Felder et al. 2002; Kallosh et al. 2002; 
Alam, Sahni & Starobinsky 2003; Frampton & Takahashi 
2003; Kallosh & Linde 2003; Sahni & Shtanov 2003). 

In a spatially flat Freidmann— Lemaitre— Robertson— 
Walker (FLRW) background, the equations of motion for a 
scalar field take the form 

3M p iH = pm + p r + p<j> , (4) 

M 2 l (2H + 3H 2 ) = -^-p4 > , (5) 

e4> + 3 He<j> + = 0, (6) 

acp 

where, 

p* = + v(0), 

1 . 2 

= 2 e< ^ -^W); «v>=?Vp4>, (7) 


which shows that, Om(z) = flom for ACDM, whereas and e = +1 and-1 corresponds to ordinary scalar and phan- 

Om(z) > fiom for quintessence (w > — 1) and Om(z) < flom tom field respectively. 
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Figure 2. This figure shows different plots for the potential m with m = 1. The upper left plot shows the evolution of energy density 
versus the redshift z. The dotted, dashed and black lines correspond to the energy density of radiation, matter and scalar field respectively. 
Initially, the energy density of scalar field (p^) is extremely sub-dominant and remains so, for most of the period of evolution. At late 
times, the field energy density catches up with the background and overtakes it. Around the present epoch, the field energy density starts 
decaying and correspondingly its equation of state starts approaching towards matter. The upper middle plot shows the evolution of 
density parameter (Q) versus the redshift z. The upper right plot shows the la (dark shaded) and 2cr (light shaded) likelihood contours 
in the — flom plane, where A i is the initial value of A. The black dot designates the best fit value. The lower left plot shows the 
evolution of q versus the redshift z. The value of q increases at redshift z < 0.4 and becomes positive at late times. The lower middle 
and lower right plots show the evolution of w(z) and Om(z) versus the redshift z. In all lower plots, dashed line represents ACDM with 
Qom = 0.315; solid (middle) line inside shaded regions show best fitted behaviour and shaded regions show lcr confidence level. We have 
used joint data (SN+Hubble+BAO+CMB), see appendix A for details. 


Let us introduce the following dimensionless quantities, 


w 


V6HM pl ’ V3HM pl ’ 


A = -M pl y, (8) 


which are used to form an autonomous system, 


where, 

H —3 — 3ea: 2 + 3i/ 2 — 
#2 = 2 
(j> _ —6x + VEey 2 X 
H<j>~ 2x 
r = VV d± 

“ ^ ’ 


(13) 

(14) 

(15) 


and prime (') denotes the derivative with respect to In a. 


/ 

X 

- x (± -JL) 

\H4> H 2 > 

(9) 


/ [3 H \ 

(10) 

V 

= -y{\l2 Xx+ ip)’ 

o! r 

= -2Q.r (2 + -^ 2 ), 

(11) 

A' 

= VQx\ 2 (1-V), 

(12) 


3.1 Quintessence field and slowing down of cosmic 
acceleration 

In this sub-section, we shall consider scalar field models 
which can facilitate the slowing down of cosmic acceleration, 
at late times. Figure |T| shows that Om(z) increases with in¬ 
creasing constant equation of state, at late times. This effect 
could correspond to dark energy decaying into dark matter 
or something else. To check this possibility, let us consider 
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Figure 3. This figure is similar to figure (2l for the potential ()l61) 
and f2om = 0.2663, where A* is the initial value of A. 

the following power law potentials, 

V(4>) = Vo4> 2m , to = 1,2 (16) 

In this case, we have slow roll regime followed by a fast roll 
and the oscillatory phase thereafter as (j> approaches the ori¬ 
gin. The average equation of state during oscillations is given 
by < w >= (to—1) / (to+ 1) and therefore, < w >= 0,1/3 for 
m = 1, 2 respectively. In this case, if we set parameters such 
that evolution of Universe corresponds to slow roll regime 
in the recent past, then acceleration would register its slow 
down around the present epoch. Keeping this in mind, we 
evolve the equations of motion numerically. Our results are 
displayed in figures [2] and [3] In figure [2J the upper right plot 
shows the la (dark shaded) and 2a (light shaded) likelihood 
contours in the A i — flom. plane, where Xi is the initial value 
of A. The best fit values of the parameters are, Xi = 1.2 and 
fi 0m = 0.2657. The lower middle and right plots show the 
evolution of w(z) and Om(z) versus the redshift 2 respec¬ 
tively. We find that the EOS of dark energy grows at late 
times which is clearly shown by Om and q plots. This corre¬ 
sponds to an increase of Om and q at redshift 2 < 0.4. We 
have used, SN+Hubble+BAO+CMB joint data for analysis 
(see appendix A). In the parametric space, shown in Ai—Uom 
plane, we find that acceleration begins to slow down around 
2 ~ 0.4 followed by deceleration at present epoch (see the 
lower plots of figures [5] and [3]). We also note from Om plot 
that Om(z) has negative curvature. The negative curvature 
of quintessence allows us to distinguish it from ACDM as 


with m = 2. The best fit values of the parameters are, A i = 1.4415 


well as from phantom dark energy irrespective of a given 
current value of the matter density. 


3.2 Slowing down of late time cosmic acceleration 
and age consideration 

The role of gravity in the Universe filled with standard mat¬ 
ter is to decelerate the expansion such that larger is the 
matter density, faster will be the expansion rate thereby 
smaller would be the age of Universe. In the standard lore, 
the only way to circumvent the age problem is provided by 
invoking a repulsive effect which could be described by a 
positive cosmological constant or by a slowly rolling scalar 
field. The repulsive effect becomes important at late times 
giving rise to large contribution to the age of Universe. It is 
therefore not surprising that more than half of the contri¬ 
bution to the age of Universe comes from the redshift in the 
interval 2 £ (0,1), see figure 0] The slowing down of accel¬ 
eration also takes place in the said interval for chosen set of 
model parameters. This should clearly decrease the age of 
Universe. Then the requirement of consistency with data can 
only be achieved by adjusting the matter density parameter 
Uo m which in a sense quantifies the attractive effect of grav¬ 
ity. Indeed, in case of 2 potential, for which slowing down 
commences around 2 ~ 0.4, the best fit value of the matter 
density parameter is, Uom = 0.2657 which is smaller than 
its ACDM value quoted by the Planck (see Table [T]). Let us 
also note that the best fit value of matter density parameter 
for the model (w(z) = — 1 + tanh l( z ~ zt ) A l ) given by Shafieloo 


but 
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Figure 4. The dimensionless T(z) = Hot(z) (t(z) is the cosmic 
age at redshift z) is plotted for three models, namely, ACDM 
(solid line), quintessence with cf) 2 potential (dashed line) and 
the cosmological model with given equation of state, w(z) = 
— 1 + tanfa [( z ~ z *) A ] (dot dashed line) which has same number of 
free parameters as the CPL ansatz but does not permit the cross¬ 
ing to phantom divide at w = — 1. The form of w(z) has been 
used (Shafieloo et al. 2009) to show slowing down of cosmic ac¬ 
celeration at £ < 0.3. The best fit value of tto m = 0.3086 for 
ACDM, is taken from Planck data result 2013, and the best fit 
values of the model parameters of 0 2 potential are found to be 
A i = 1.2 and f2om = 0.2657, (where A i is the initial value of 
A) whereas the best fit values of parameters for the ansatz are, 
z t = 0.008, A = 12.8 and Q 0 m = 0.255 (Shafieloo et al. 2009). 
The horizontal band represents the range of the age of Universe 
based upon observations on globular clusters (Krauss &; Chaboyer 
2003; Frieman, Turner &; Huterer 2008). The black dot designates 
the numerical value of T(z) at z = 1 ( T(z = 1) = 0.5498). 


Table 1. Age summary. 


Data/Model 

^Om 

Ho 

(Km/s/Mpc) 

to 

ipy) 

Planck 

0.3086 

67.77 

13.7965 

ACDM 

0.3086 

67.77 

13.8567 

cf) 2, potential 

0.2657 

67.77 

13.4825 

w _ l+tanh[(z-z t )/\\ 

0.255 

67.77 

14.4039 


et. al (2009), is smaller than its counter part obtained in 0 2 
model giving rise to some what larger value of the age (see 
figure |4j). 

The cosmic age t(z) is defined as the cosmic time, Uni¬ 
verse spends from a given 2 to the present epoch as 



dz' 

(1 + z')H(z') ’ 


(17) 


where H(z) is the Hubble parameter expressed in terms of 
redshift z. For spatially flat Universe, Hubble parameter can 
be written as 

H(z) = tfo[fWl + z) 3 + (l-fl 0m ) 

exp ( 3 / [ 1 + w( U] 1T7)] 7 ’ (18) 


where w = —1 for ACDM, w(z) = -ri— tJt — for quintessence 

with <j> 2 potential and w(z) = — 1 + tanh l( z ~ z t') A l f or the cos _ 
mological model suggested by Shafieloo et al. (2009). 

It is helpful to use a dimensionless cosmic age T(z) = 
Hot(z), where t(z) is given by equation (1 1 71) such that 
T(z) = Hot(z) 0 as z —>• 0 implying (a) t(0) —> 0, (b) 
1(0) <C t(z > 1) as shown in figure U 


3.3 Tracker 

Scalar field models with tracker solutions are of great impor¬ 
tance in cosmology. A tracker solution is such that, during 
radiation and matter era, the field mimics the background- 
scaling regime and at late times it exits to late time cosmic 
acceleration. Once the model parameters are fixed, the exit 
to dark energy behaviour remains independent for a wide 
variety of initial conditions. Here we shall consider the fol¬ 
lowing potential (Sahni & Wang 2000): 

V{<f>) = Vo [cosh (&c/>/M p ) - l] p , 0 < p < 1/2 (19) 

The above potential has asymptotic form 

V{4>) = » 1 , (20) 

and 

^ (jfy P , P&W/Mp « 1 , ( 21 ) 

where pa = a. The power law behaviour of m near the 
origin leads to oscillations of <f) when it approaches the origin. 
To have viable thermal history of the Universe, we need 
to have = constants 0.01 (Ade et al. 2013) 

during the radiation dominated era, which implies a > 20 
(here wb represents background equation of state). In this 
case, the scalar field behaves like background matter, for 
most of the history of Universe and, only at late times, it 
exits to dark energy which demands for tracker solution. 

The potential under consideration has peculiarity near 
the origin for generic values of p. We have noticed that it 
is problematic to obtain oscillatory behaviour in the convex 
core of m, using the autonomous system of equations. In 
what follows, we shall use the dimensionless form of evo¬ 
lution equations which allows us to control the dynamics 
of field, beginning from scaling regime to oscillations about 
<j> = 0. The equations of motion have the following form, 

^ 4 ! = smi [p * +pm+pr] (22) 


(}i + 3—<j> -f- Vfi^) — 0. 
a 


(23) 


In order to investigate the dynamics described by equations 
(1221) and (12311 . it would be convenient to cast them as a 
system of first order equations 


Ylm 


Yo 


h(Y u Y 2 ) 


(24) 


Yi = —3Y 2 - 


1 

h(Y 1 ,Y 2 ) 


rdV(Yi) 
L dYi 


(25) 
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ln(l+z) 


Figure 5. This figure shows the evolution of potential (1191 ) and w versus field <j> and redshift z respectively. Figure shows that at late 
times the equation of state oscillates between zero and —1 such that the system spends most of the time around w= —1. We have taken 
n 0 m — 0-3, p = 0.1, a = 200 (a = pa = 20) and Vq = 2.2 M^Hq. 




Figure 6. This figure shows different plots for the potential (|19|) . The left and middle plots show the tracking behaviour of scalar field 
energy density and evolution of density parameter versus the redshift z. The dotted, dashed and black lines correspond to the radiation, 
matter and scalar field respectively. In both the plots the parameters are same as in figure [5] The right plot shows the evolution of 
average equation of state <w> versus cosmic time which agrees with <w> = our plot corresponds to p = 1, Qom = 0.319, a = 13, 
Vq — 2.2 MpH q and p = 0.096, Qom — 0.319, a = 37.59, Vq = 2.2 M^Hq. The vertical dotted line designates the present epoch. 



Figure 7. This figure shows the la (dark shaded) and 2 a (light shaded) likelihood contours for the potential dj 19D with joint data 
(SN+Hubble+BAO+CMB); black dots designate the best fit value of the parameters which are found to be 9om — 0.319, p = 0.096 and 
d = 37.59. 
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Figure 8. This figure shows the evolution of w(z) and Om(z) versus the redshift z for the potential (119D . Equation of state exhibits 
oscillating behaviour which has also been manifested in Om(z) plot. In both the plots, dashed line represents ACDM with Do m = 0.319; 
solid (middle) lines inside shaded regions show best fitted behaviour and shaded regions show lcr confidence level. We have used joint 
data (SN+Hubble+BAO+CMB) in our analysis. 


where 


= A. = v = ^V 

M p ’ M P H 0 ’ 


(26) 


and prime denotes the derivative with respect to the variable 
N = ln(a). The function h(Y\, Y 2 ) is given as: 


h(Y!,Y 2 ) 


Yj V(Xi) 

6 3 


T Do mC ^ a + Uo r e ^ a 


(27) 


where Do m and Do are the present energy density param¬ 
eters of matter and radiation respectively. We have numer¬ 
ically solved the evolution equations m & 425)1; Our re¬ 
sults are shown in figures El ED CZ) an d ID Figure [5] shows 
the behaviour of the potential m versus field <j> and the 
equation of state parameter w versus the redshift a. As the 
field evolves from the steep region towards the origin, p$ un¬ 
dershoots the background and field freezes for a while due 
to Hubble damping (see the left and the middle plots of fig¬ 
ure [SJ. The field then approximately mimics the background 
before approaching the convex core of the potential, where 
oscillations set in. It is clear from the right plot of figure [5] 
that field oscillates most of the time near w= — 1. The right 
plot of figure [6] shows average value of w versus the cosmic 
time for p = 1 & p = 0.096 which agrees with the analytical 
result, <w>= p — 1/p + 1 at the attractor point. 

In figure [3 all plots show the lcr (dark shaded) and 
2cr (light shaded) likelihood contours. Figure [8] shows the 
oscillating behaviour of equation of state and the corre¬ 
sponding behaviour of Om inside lcr confidence level. The 
joint data (SN+Hubble+BAO+CMB) was used for carrying 
out the observational analysis (see appendix A). We should 
note that the best fit value for the parameter p is different 
from p = 0 that would correspond to the case of cosmo¬ 
logical constant. The latter is reinforced by figure [8] which 
clearly shows, the model under consideration differs from 
ACDM within lcr confidence level. In this case, we should 
further check for Xred- We find that the numerical value of 
Xred. = 0.849713 corresponding to the best fit value of p, is 
much smaller than one (see Table [2|. Thus, we can not con¬ 
clude that the underlying model is preferred over ACDM. 


Table 2. Comparison between Xred. values for the potential cnni 
and ACDM. To obtain X^. ed for the potential d 19 D , we vary p and 
put the best fit values of the remaining parameters, Do m = 0.319 
and a = 37.59. The bold values correspond to the best fit value 
of p and Dom for the potential on and ACDM respectively. 


Potential (19) 

V xled 

ACDM 

Do m Xred 

0.01 

0.951237 

0.25 

0.956061 

0.02 

0.973990 

0.26 

0.948286 

0.03 

0.987693 

0.27 

0.943602 

0.04 

0.992347 

0.28 

0.94174 

0.05 

0.985540 

0.29 

0.942457 

0.06 

0.972094 

0.30 

0.945545 

0.07 

0.952011 

0.31 

0.950822 

0.08 

0.888724 

0.32 

0.958119 

0.09 

0.855366 

0.33 

0.967293 

0.096 

0.849713 

0.34 

0.978213 

0.10 

0.851936 

0.35 

0.990760 


In this case, we could not find slowing down effect, for any 
values of the model parameters. 


3.4 Phantom Held with linear potential 

Phantom dark energy with equation of state w < — 1 can be 
achieved by introducing a negative kinetic energy term in the 
action of the scalar field. By putting e = — 1 in equations 
© and 0, we get equation of motion, energy density and 
pressure of phantom field (Singh, Sami & Dadhich 2003). In 
what follows, we shall examine the dynamic of a phantom 
field. In this case, the Hubble parameter for spatially flat 
Universe can be written as 

H\z) = H 0 2 [n 0 r(l + 2) 4 + Dom(l + +) 3 

+D 0 «/,exp^3 [l + w^z')] ^ ^ ], (28) 

where Do, Dom and Do</> are the present energy density 
parameters of radiation, matter and field respectively, and 
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Figure 9. This figure shows different plots for the linear potential. The upper left plot shows the evolution of energy density versus 
the redshift z. The dotted, dashed and black lines correspond to the energy density of radiation, matter and phantom field respectively. 
Initially, the energy density of phantom field is extremely sub-dominant and remains to be so, for most of the period of evolution. At 
late times, the field energy density catches up with the background, overtakes it, and starts growing (w < —1) and derives the current 
accelerated expansion of the Universe. The upper right plot shows the la (dark shaded) and 2 a (light shaded) likelihood contours in 
the A i — f2om plane (the black dot designates the best fit value of A i and fiom, where A i is the initial value of A). The lower left and 
right plots show the evolution of w(z) and Ora(z) versus the redshift z respectively. Om(z) has positive curvature. In both lower plots, 
the horizontal dashed line represents ACDM with f2o m =0.3, solid (middle) line inside shaded regions show best fitted behaviour and 
shaded regions show lcr confidence level. We have used joint data (SN+Hubble-bBAO+CMB) in carrying out the analysis. 


Ho designates the value of Hubble parameter at the present 
epoch. We will be interested in the phantom dynamics with 
a linear potential, 


Vtf) = Vo<t>, (29) 


Numerically integrating the equations of motion, we find the 
field energy density, equation of state and Om for the said 
potential. The results are shown in figure [HI The upper left 
plot shows the evolution of energy density versus redshift 
s whereas the upper right plot shows the lcr (dark shaded) 
and 2cr (light shaded) likelihood contours in the Xi — fiom 
plane. The best fit values of the parameters are found to 
be A i = 0.1959 and flom = 0.2826. The lower plots show 
the evolution of w(z) and Om{z). As seen in the figure, Om 
has positive curvature for phantom field model which distin¬ 
guishes phantom field model from zero-curvature (ACDM), 
for any current value of the matter density. The joint data 
i.e. SN+Hubble+BAO+CMB was used for analysis (see ap¬ 
pendix A). 


4 NON MINIMALLY COUPLED SCALAR 
FIELD MODEL 

In this section, we revisit a non minimally coupled scalar 
field model that allows to obtain a transient phantom dark 
energy. Various features of non minimally coupled scalar 
field system have been investigated (Dent et al. 2013; Ka- 
menshchik et al. 2013; Nozari & Rashidi 2013; Aref’eva et 
al. 2014; Kamenshchik et al. 2014; Luo, Wu & Yu 2014; 
Pozdeeva & Vernov 2014; Skugoreva, Toporensky & Vernov 
2014; Skugoreva, Saridakis & Toporensky 2014). The action 
having non minimal coupling with scalar field is given as 
(Sami et al. 2012): 

5 = \J V=gd*x [| - (<^<M„ + £RB(<I>) + 2V(</>))] +S M , 

(30) 

where k=8itG, £ is the dimensionless coupling constant and 
Sm is the matter action. The equations of motion which are 
obtained by varying the action HI have the form (Sami et 
al. 2012): 

H 2 = | + V(4>) + 3 + H 2 B{4,)) + p) , 


( 31 ) 
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Figure 10. This figure corresponds to the case of non minimal coupling with N = 2, n = 7, £ = 0.2 and Qo m — 0.316. The left and 
right plots show the evolution of (1 — 6£13(</>)) and w versus the redshift z respectively, where G e ff = 6/87r(l — .£?((/>)). The black dots 

( point on the vertical dashed line) designate the present epoch which occurs in the regime of G e f / > 0, the effective Newtonian constant 
changes sign thereafter in future. The point on the horizontal dashed line is the epoch where G e ff turns negative. 



Figure 11. This figure represents the non minimally coupled scalar field DE model with IV = 2, n = 7; left plot shows the la (dark 
shaded) and 2cr (light shaded) likelihood contours in the £ — Qom plane. The black dot corresponds to the best fit value. Evolution of 
w(,z) is shown in the middle plot. The right plot shows the evolution of Om(z) versus the redshift z. In middle and right plots, solid 
(middle) line inside shaded regions show best fitted behaviour and shaded regions show la confidence level; the horizontal dashed line 
represents ACDM with Qom = 0.315. The SN+Hubble+BAO+CMB data has been used jointly in carrying out the analysis. 


R = k( - 0 2 + + 3£(3 H<j>B’(<j>) + jB{<j>) 

+ (</,)) + p(l-3w)), (32) 

j> + 3H<j>+ltRB'(<t>) + V'(<t>) = 0, (33) 

where R = 6(2 H 2 + H), is the Ricci Scalar and p is the 
energy density of the matter. 

In this case (Sami et al. 2012), we have a varying ef¬ 
fective Newtonian gravitational constant G e // which is a 
function of scalar field <j>. We found stationary points and 
their stability by considering B((f)) = (j> and V{(j>) = Voijf 1 . 
The de Sitter solution of interest to us, exists for n < 2 N 
for which G e // = 6/87r(l — 6 £B((p)) = 3(2 N — n)/8nN is 
positive. In order to check the stability of the de Sitter so¬ 
lution, we consider small perturbations, and get the system 


of equations (see appendix B) 



The numerical analysis of the system exhibits that the de 
Sitter solution is stable when 2 N + 1 ^ n and for positive 
values of the coupling, £ > 0 (we assume Vo = 1). Thus 
we display few cases in Table Q5] which show how the na¬ 
ture of the fixed points depends upon numerical values of 
As demonstrated by Sami et al. (2012), the de Sitter so¬ 
lution in this case occurs in the region of negative effective 
gravitational constant, thereby leading to a ghost dominated 
Universe in future and a transient quintessence (phantom) 
phase with G e // > 0 around the present epoch. Figure fTTil 
shows that before going to de Sitter point with G e // < 0, 
the equation of state passes through a phantom phase. In 
order to obtain phantom (w < —1) phase consistent with 
observation and G e // > 0 at present epoch, we adjust the 
model parameters as shown in figure [TOl 
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Table 3. Nature of the fixed points (Sami et al. 2012): 


N n £ nature of fixed points 


2 

5 

0 < £ < 0.1333 

saddle 



0.1333 < £ < 0.2068 

attractive focus 



0.2069 < £ < 1 

attractive node 

2 

7 

0 < £ < 0.0952 

saddle 



0.0952 < £ < 0.3999 

attractive focus 



0.4000 < £ < 1 

attractive node 

2 

9 

0 < £ < 0.0740 

saddle 



0.0740 < £ < 1 

attractive focus 

4 

9 

0 < £ < 0.0009 

saddle 



0.0009 < £ < 0.0014 

attractive focus 



0.0015 < £ < 1 

attractive node 


Let us now apply Om diagnostic to the case under con¬ 
sideration and examine the curvature (slope) of Om. To this 
effect, we investigate the evolution equations numerically in 
case of IV = 2 and n = 7 for a viable range of parameters 
£ and fi 0 m (0.24 < fl 0 m < 0.36 and 0.12 < £ < 0.61); the 
range of f; is dictated by stability considerations. Our results 
are displayed in figure fill in which the left plot shows the 
la (dark shaded) and 2a (light shaded) likelihood contours 
in the £ — flom plane. The best fit values of the parameters 
are, £ = 0.4574 and flom = 0.3393. The middle plot shows 
the evolution of w(a) versus redshift 2 whereas the right plot 
displays the evolution of Om(z) versus redshift 2 . We find 
that Om has positive curvature (slope) for phantom equa¬ 
tion of state which is a generic feature of dark energy models 
with w < — 1. The positive curvature of phantom distin¬ 
guishes non minimally coupled scalar field model from zero 
curvature ACDM model, for any value of the matter density 
as shown in figure 1111 We used SN+Hubble+BAO+CMB 
joint data for carrying out the observational analysis, see 
appendix A for details. 


5 CONCLUSIONS 

In this paper, we have investigated scalar field models (in¬ 
cluding the phantom case) using the Om diagnostic. We 
have specifically focused on models with power law poten¬ 
tials that lead to the slowing down of late time cosmic 
acceleration. In case of quintessence with quadratic and 
quartic potentials, we have demonstrated that the slowing 
down phenomenon takes place for 2 < 0.4. Figures [2] and 
[3] show that at late times, the field energy density starts 
decreasing (w > —1) and exhibits matter/radiation like be¬ 
haviour corresponding to the deceleration of expansion. Con¬ 
sequently, the equation of state of dark energy, Om(z) and 
the deceleration parameter q grow for redshift in the inter¬ 
val 2 £ (0, 0.4). This signifies that cosmic acceleration might 
have already peaked and that we are presently observing its 
slowing down. We find that quintessence models have neg¬ 
ative curvature that differentiates these models from zero 
curvature (ACDM), for any given current value of the mat¬ 
ter density as shown by the Om plots in figures [2] and [3j 
The best fit values of the model parameters for (f ) 2 and 0 4 


potentials are found to be A; = 1.2, flo m = 0.2657 and 
A i = 1.4415, flom = 0.2663 respectively. In these models, 
the best fit value of flom is always less than that for its 
counter part i.e. ACDM, which compensates for the effect 
of intermediate slowing down. We have also investigated a 
model with a cosh potential GHJ. which has the tracking 
property. For small values of 4>, the said potential mimics a 
power law behaviour and gives rise to oscillations of 4> near 
the origin. The tracking behaviour of the scalar field energy 
density, the evolution of density parameter, the oscillating 
behaviour of the equation of state at late times and the cor¬ 
responding behaviour of Om ( within la confidence level ) 
are shown in figures [6] and [8] The field oscillations in the 
convex core around the origin are reflected in the behaviour 
of w (see figure [SJ), such that the average equation of state 
parameter is given by <w>= p — 1/p + 1 (see figure [6]). 
The best fit values of the model parameters are, p = 0.096, 
a = 37.59, fiom = 0.319. In this case, the ACDM is clearly 
outside 2a confidence level as shown in figure [3 This is also 
clear from the figure [8] and it deserves a comment. In order 
to draw a final conclusion, we looked for Xred, and found that 
Xred is much smaller than one. Thus we can not claim that 
the model under consideration is better than ACDM. 

Next, we applied the Om diagnostic to phantom dark 
energy. In this case, we considered the phantom field with 
a linear potential. In figure [9] we have shown that Om has 
positive curvature that distinguishes the phantom dark en¬ 
ergy from the zero curvature ACDM, for any current value 
of the matter density. The best fit values of the model pa¬ 
rameters are, A i = 0.1959, flom = 0.2826. We also examined 
a non minimally coupled scalar field model, which has a 
transient phantom behaviour. In this case again, Om(z) has 
positive curvature as shown in the figure 1111 The best fit 
values of the model parameters are found to be £ = 0.4574, 
Dom = 0.3393. It signifies that the positive curvature of Om 
is a generic feature of phantom dark energy, be it transient 
or otherwise. 

We conclude that given the present data, the Om di¬ 
agnostic can clearly distinguish between scalar field models 
and ACDM and that in case of quadratic and quartic po¬ 
tentials, there exists a specific region in the parameter space 
which could allow for the slowing down of late time cosmic 
acceleration. 
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APPENDIX A: OBSERVATIONAL DATA 
ANALYSIS 

We put constraints on the model parameters using recent 
observational data, namely Type la Supernovae, BAO, CMB 
and data of Hubble parameter. The total \ 2 f° r joint data 
is defined as 

Xtot = XSN + XbAO + XHub + XCMB ; (35) 







11 


Table 4. Values of d A V *) — f or distinct values of zrao- 
Dy(Z BAO ) ^ 


ZbAO 

0.106 

0.2 

0.35 

0.44 

0.6 

0.73 

d A(z*) 

D v (Zbao) 

30.95 ± 1.46 

17.55 ±0.60 

10.11 ±0.37 

8.44 ±0.67 

6.69 ± 0.33 

5.45 ±0.31 


where the \i f° r each data set is evaluated as follows: First, 
we consider the Type la supernova observation which is one 
of the direct probes for the cosmological expansion. We use 
Union2.1 compilation data (Suzuki et al. 2012) of 580 data 
points. For this case, one measures the apparent luminosity 
of the supernova explosion from the photon flux received. In 
the present context, one of the most relevant cosmological 
quantity is luminosity distance Dl(z) defined as, 

-(l + *>/‘gg- (36) 

Cosmologists often use the distance modulus /x defined as 
p = m — M = 51ogDi, + /xo, where m and M are the 
apparent and absolute magnitudes of the Supernovae and 

po = 5 log ^ ^p C ^ + 25 is a nuisance parameter which is 

marginalized. The corresponding \ 2 is written as 

.2 ( .. a\ _ \_Pth(Zi, Ho, 9) — p,obs (Zi)] f ^ 

Xsn(M o, V) - 2^- ( 7 .\2 -> l 37 ! 


where p, 0 bs, pth and represents the observed, theoreti¬ 
cal distance modulus and uncertainty in the distance mod¬ 
ulus respectively; 9 represents any parameter of the partic¬ 
ular model. Eventually, marginalizing /xo following Lazkoz, 
Nesseris & Perivolaropoulos (2005), we get 


Xsn (6) 


m 


B(9) 2 

w 


(38) 


where, 


m 

B{9) 

C{6) 


580 

E 

i=1 
580 

E 

i= 1 
580 


E 


[fJ'th (Zi , /Xo — 0,6) /Xo6s(*i)] 

Vli{Zi) 2 

pthi^Zi, /X0 — 0, 9) pobsi.Zi) 

’ 



(39) 

(40) 

(41) 


Next, we use BAO data of D d fz Z BA ) (Eisenstein et al. 2005; 
Percival et al. 2010; Beutler et al. 2011; Blake et al. 2011; 
Jarosik et al. 2011; Giostri et al. 2012), where 2 * « 1091 
is the decoupling time, d A (z) = h(z') co-moving 

angular-diameter distance and Dv(z) = ( '^d A (z ) 2 H z ^ j ? is 
the dilation scale. Data required for this analysis is shown 
in Table U 

The Xbao is defined as (Giostri et al. 2012), 


(42) 


where, 


/ 


X = 


V 


d A G*) 

D v (0.106) 

d A<.~*) _ 

Dy (0.2) 
d A ( z +) 

Dy (0.35) 

d A(-*) 

Dy (0.44) 
d A(z*) 

Dy (0.6) 
d A(z*) 

Dy( 0.73) 


- 30.95 
17.55 
10.11 

- 8.44 
6.69 

- 5.45 


\ 


/ 


(43) 


and C 1 is the inverse covariance matrix given by Giostri et 
al. (2012). 

We then use the observational data on Hubble param¬ 
eter as recently compiled by Farooq & Ratra (2013) in the 
redshift range 0.07 < 2 < 2.3. The sample contains 28 ob¬ 
servational data points of H(z). These values are given in 
Table [5] To complete the data set, we take the latest and 
most precise measurement of the Hubble constant Ho from 
PLANCK 2013 results (Ade et al. 2013). To apply the data 
of Hubble parameter on our models, we work with the nor¬ 
malized Hubble parameter, h = H/Hq. 

The x 2 for the normalized Hubble parameter is defined 
as, 


29 

XHub(#) = E 


i=l 


[hth (Zi,9) hobs ( V )] 
CTh(Zi) 2 


(44) 


where, h 0 o s and hth are the observed and theoretical values 
of the normalized Hubble parameter respectively. 

Also, 




(j Hq \ 

V H ^ Ho ) 


h, 


(45) 


where <jh and <jh 0 is the error in H and Ho respectively. 

Finally, we apply CMB shift parameter R = 
HoV^mo fg 089 b(z') ■ Tlle corresponding Xcmb can be writ¬ 
ten as, 


Xcmb(#) 


(R(9) - Rg) 2 


(46) 


where, Ro = 1.725 ± 0.018 ( Komatsu et al. 2011). 


APPENDIX B 


For de Sitter solution, substituting H = <j> = <j> = p = 0 and 
R = 12Ho 2 with B(rj)) = cf> N and V(4>) = Vo<f> n in equations 
(ED, ED and ED, we obtain 


H 0 \1-6^o N ) = 2Vo4>o n , 

QH 0 2 (,N(j)o N ~ 1 ± Vonfo”- 1 = 0 , 

which gives 

o _ V 0 n<j> 0 n - N N _ n 

0 6 £N ’ 6£(n - 2 N )' 


(47) 


(48) 


The de Sitter solution exists provided that n < 2N, and 
G e ff = 6/87t(1 — 6£.B(<)>)) = 3(2 N — n)/8nN is positive 


xIao =X t C~ 1 X, 
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Table 5. H(z) measurements (in unit [kms x Mpc *■]) and their 
errors (Farooq & Ratra 2013). 


z 

H(z) 

°H 

Reference 

0.070 

69 

19.6 

Zhang et al. 2012 

0.100 

69 

12 

Simon et al. 2005 

0.120 

68.6 

26.2 

Zhang et al. 2012 

0.170 

83 

8 

Simon et al. 2005 

0.179 

75 

4 

Moresco et al. 2012 

0.199 

75 

5 

Moresco et ah 2012 

0.200 

72.9 

29.6 

Zhang et al. 2012 

0.270 

77 

14 

Simon et al. 2005 

0.280 

88.8 

36.6 

Zhang et al. 2012 

0.350 

76.3 

5.6 

Chuang & Wang 2012 

0.352 

83 

14 

Moresco et ah 2012 

0.400 

95 

17 

Simon et al. 2005 

0.440 

82.6 

7.8 

Blake et al. 2012 

0.480 

97 

62 

Stern et al. 2010 

0.593 

104 

13 

Moresco et ah 2012 

0.600 

87.9 

6.1 

Blake et al. 2012 

0.680 

92 

8 

Moresco et al. 2012 

0.730 

97.3 

7.0 

Blake et al. 2012 

0.781 

105 

12 

Moresco et ah 2012 

0.875 

125 

17 

Moresco et ah 2012 

0.880 

90 

40 

Stern et al. 2010 

0.900 

117 

23 

Simon et al. 2005 

1.037 

154 

20 

Moresco et ah 2012 

1.300 

168 

17 

Simon et al. 2005 

1.430 

177 

18 

Simon et al. 2005 

1.530 

140 

14 

Simon et al. 2005 

1.750 

202 

40 

Simon et al. 2005 

2.300 

224 

8 

Busca et al. 2013 


only if n < 2N. The system of equations m, ® and (|33l> 
for p = 0 then reduces to, 

H(1 - 6 £-B + 9£ 2 B' 2 ) = 30 2 (fB" - 1) 

-3£B'{4H<j> + V' + 6{B'H 2 ), 
4>(1 - 6 £-B + 9£ 2 B' 2 ) = -(3H<j>+V'){l-6£B) 

-£B'(-3<j> 2 + 12V + 9£(3 H<j>B' 
+<j> 2 B")). (49) 

To check the stability of the de Sitter, we consider small 
perturbations p & v around this background: H = Ho + p 
and cj> = 4> o + ^ in the dynamical system flSD> which gives us 
the evolution equations for perturbations 

p(36£ 2 N 2 H 0 cl>o 2N - 2 ) 

+p( 1 - 6^o N + 9li 2 N 2 cf> 0 2N - 2 ) 

+is(3£V 0 Nn(N + n- 2)0 O JV+71_3 
+36£ 2 N 2 (N - 1)H 0 2 <(>o 2JV " 3 ) 

+z>(12£AhH 0 <(>o JV “ 1 ) = 0 , (50) 

v(Von(n — l)0 o n_2 (l — &£,<j>o N ) 

—6£,VoNnrj>o N+n ~ 2 + 12^V 0 N(N + n- 1 )</> 0 JV+n “ 2 ) 
+!>(3JT 0 (1 - 6£<t> 0 N ) + 27fN 2 H 0( p 0 2N - 2 ) 

+/>(!- 6 ^ 0 ^ + H 2 N 2 <i>o 2N - 2 ) = 0. (51) 


Equations m and m can be put in a simple form by 
introducing the following notations 

+i = 36ti 2 N 2 H 0 (j> 0 2N - 2 , 

= 1 -&ti<t>o N +VZ 2 N 2 <t>o 2N - 2 , 

^3 = 3£V 0 Nn{N + n- 2 )<t> 0 N+n ~ 3 

+36(, 2 N 2 (N — 1)Hq 2 <j>o 2N ~ 3 ^ 

A 4 = U^NHofo"- 1 , 

A 5 = Von(n — l)</> 0 n-2 (l — 6£(f>o N ) — 6^VoNnrj)o N+rL ~ 2 
+12 £V 0 N(N + n- 1 )<(>o JV+rl “ 2 , 

A 6 = 3JTo(l-6^o JV ) + 27e 2 A r2 -H r o<(>o 2iV_2 , 

A 7 = 1-6^o N +9{ 2 N 2 0 o 2N ~ 2 (52) 

Next, by using, a = p, /3 = v, 7 = v, the system of equations 
(|50ll and m acquires a simple form, 

A\ol + A 2 P + As/3 + A 4 ^ = 0, (53) 

A 5 fi + +67 + A 7 v = 0. (54) 

Taking derivative of a, (3 and 7 with respect to time we get 
the system of equations, 
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